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ABSTRACT 

Cooperative effects arising from the contributions of a macroscopic assem- 
blage of interacting sites will greatly magnify the influence of weak-field 
interaction terms. As a result, significant differential adsorption of 
optical isomers onto an achiral surface is predicted to occur if the surface 
is continuous and sufficiently large. 

I. Introduction 

The possibility that achiral surfaces, such as clay surfaces, may still 
differentiate between optical isomers has intrigued researchers and has been 
the motivation for numerous investigations [1,2, 3, 4]. The results to date 
have been ambiguous. A theoretical framework, which could be used to decide 
whether or not a set of experiments might be successful or even delineate 
the minimal conditions necessary for success, has been lacking. It is the 
purpose of this communication to provide a possible theoretical model for 
the resolution of enantiomers by achiral surfaces. 

^National Academy of Sciences Research Associate. 
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The question of chirality has been subjected to reexamination recently 
in two very different approaches. The first adopts Hund's explanation [5] 
but appends weak-field interactions to it. The second rejects Hund's explana' 
tion and replaces it with a possibly more fundamental one [6]. Weak-field 
interactions are not included in Pfeifer's [6] treatment, but such inclusion 
is thought to be straightforward. The importance of the weak-field inter- 
actions are that they break the degeneracy between optical isomers to a 
miniscule extent, generally. We will argue that there are cases where the 
effects will not be insignificant. Thos the opening remarks in this para- 
graph are important in establishing that effects arising from weak-field 
interactions are not dependent on the particular model of chirality. The 
extent to which the degeneracy between optical isomers is split may be model- 
dependent, however. The discussion of the weak-field interaction begins with 
the formal expression obtained by Bouchlat and Bouchiat [7]. 


II. Weak-Field Interaction Term s 

Nonvanishing contribution to the parity nonconserving energy arises in 
second order when spin-orbit coupling is included. Approximating the ground 
state and excited wavefunctions of a molecular system by single determinants 
of molecular orbitals yields, following Rein et al. [8], the pnc energy 
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The (x) terms are spin functions while f denotes filled orbitals and e 
denotes empty ones. Summing over the spin magnetic quantum numbers In 
equation (1) reduces the expression to the spin-independent form (from Rein 
et al. [8]) . 


6w ■ 21 ■ 


“'“7 — + c ’ c - 


( 2 ) 


Here £ and £ are the appropriate operators obtained from H pnc and H 80 , respec- 
tively, upon suppression of the spin operators. The molecular orbitals are 
assumed to be expressible as LCAOs where 

♦f ’ i k(« C * M»*HU <3) 

The are atomic orbitals associated with the 1th atomic center and the 

Cf are the appropriate coefficients of the in the f molecular 

orbital. Again following Rein et al. [8] , overlap is neglected and each 
matrix element in equation (2) is restricted to a single center. Then, with 
use of equation (3) , equation (2) yields 
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As Rein, et al. noted [8] , 6w - 0 if 1 - r. The importance of the expres- 
sion in equation (4) is Its generality. It Is net restricted to simple 
molecules but can be applied equally well to extensive multicenter inter- 
acting systems such as crystals. For such a system, the atomic orbitals from 
which the molecular orbitals are constructed are dispersed throughout and 
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include the entire crystal. Thus equation (4) , in this context, leads to the 
implication that cooperative effects are operating in a giant molecule with 
many sites. As a coneequence, all sites (assuming that each site experiences 
some degree of asymmetry in its environment) contribute to the energy 6ui 
in such a fashion that there can be no reduction of the double summation 
over those sites. The operators N and A are preferentially weighted toward 
those regions close to "heavy" nuclei. Accordingly, if all sites are identical 
and each contains one "heavy" atom (such as silicon) and there are N sites, 
then, for those molecular orbitals whose composition includes orbitals 
associated with the ith "heavy" nuclear center, the coeffients of these 
atomic orbitals are 


S.moI ' N ' 



Furthermore, the sums in equation (4) for these molecular orbitals run over 
all N "heavy" centers. The number of molecular orbitals that contain atomic 
orbitals associated with the "heavy" nuclear centers is 0(N). Therefore, 
the weighted contribution from all the sites to a <•••*<•••> term in 
equation (4) is 0(1) . 

The unperturbed energy for the system, relating it to an energy, e, 
associated with a hypothetical isolated site, is set equal to 
u « Ne. 

Similarly, if we denote by 6c the value of the average weak-field interaction 
term that can be assigned to the hypothetical isolated site, examination of 
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equation (A) suggests that the second-order perturbation term equals 
6w « N 2 fie. 

Any proportionality constat hss *u*eri absorbed in the fie. This relationship 
is supported upon noting that the sum over f, e orbitals in equation (A) 
Involves 0(N 2 ) terms. Substitution of the denominator in equation (A) by 
< Aw > e j, where < Aw > e j is the weighted difference of a representative 
subset of (id, - w ) terms from the f and e bands, yields 

X 6 
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Here the (£*, e*) sum is over an appropriate subset of f and e orbitals; 
these have been labeled f* and e*. Thus there appears to be an enhancement 
effect of 0(N) in a system comprised of N coupled sites. Consequently, this 
imposes a constraint on the second-order perturbation expression since 
equation (A) can only be assumed valid if w » fiw. (In terms of our approxi- 
mate expressions, e » Nfie.) The following section will focus on the applica- 
tion to very large systems by using the techniques of statistical mechanics. 


III. Application of the Weak-Fleld Interaction Perturbation 

The magnitude of the weak-field terms for an isolated site is very 
small; estimates range from 10“ 21 au to 10“ 18 au [5,8], Therefore, although 
it has been argued [9] that such infinitesimal differences cannot affect 
chemical processes, cooperative effects will profoundly alter this conclu- 
sion, We adopt a picture in which adsorption occurs from a solution 
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containing a racemic mixture (the optical iaomers are hereafter noted as A and 
B) onto a continuous surface possessing N achiral sites, each of which has 
a "heavy" atom. It is assumed that all sites are occupied with either an 
A or a B isomer. The spin-orbit terms in equation (4) include all N sites. 

Two cases are considered. In the first case we assume that the adsorbed 
isomer introduces asymmetry into the environment of the "heavy" atom and it 
is also strongly coupled with the surface. In the second case, which is a 
limit of the first, we only assume that an asymmetry is introduced into the 
"heavy" atom environment. The first case will be treated in some detail 
since the results can be easily modified to yield the second-case limit. 

The energy difference between adsorption of the A Isomer over the B form is 
given by 


AE * N(fie A - fie B + fie? - fie? ) 
site site isomer isomer 
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It is readily shown, if we assume that the Boltzmann distribution is appli- 
cable, that the partition function for this system equals 
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Here, @ - 1/kT where k is the Boltzmann constant and N is constrained by the 
requirement that 

N - N A + N b . (8) 

B B 

The terms AE, $ e 8;lte > ® e i soiner » N are constants. Only and Ng can vary; 

however, from the relationship in equation (8) the variation of is 

5N a - - 5N b . (9) 
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The most probable distribution of A and B Is found by setting the variation 

£> 

of Q, 

6Q - 0 

and solving the resulting equation. This distribution, In terms of and 
Ng» Is found to equal 

N b /N a - exp(UAE) . (10) 

Suppose there are M Isolated surfaces, each possessing N sites, which 
are simultaneously exposed to the same racemic solution. The most probable 
distribution is obtained by using the same procedure outlined previously. The 
result is identical to the relationship in equation (10), i.e., 

MN b /MN a - N b /N a - exp(SAE) . - (11) 


Thus the most probable distribution depends on the number of sites per surface. 
It depends on the total surface only If that total represents a single con- 
tinuous surface. This result might be crucial in determining whether or not 
an experiment can yield detectable differential adsorption. 

As an illustration, we assume that 


6e 

5e 


A 

site 

A 

isomer 


10 -23 au, 

a i0”20 aUt 


A 

The value for ^ e sit;e is chosen on the assumption that the estimate by Rein 
et al. [8] is representative for the situation in which asymmetry is 
introduced at a distance from the "heavy" atom. The value $ e isomer i8 P ur " 
posefully chosen two orders smaller than the upper-bound estimate (10~ 18 au) 
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for comparable ays terns since this upper bound Is In question [8] . Then the 
number of sites a single surface (hereafter crystal and surface will be used 

t 

interchangeably) must possess in order that there be a one-percent differentia- 
tion in the adsorption of B isomer with respect to A isomer can be determined 
by using equation (10). We desire that the ratio of the number of adsorbed 
A Isomers to the number of adsorbed B isomers equal 

N b /N a - 1.020. 

Solving equation (10) for BAE yields 
BAE “ 0.020. 

Since 

AE « 2N x 10" 20 au. 


then at 300°K the number of sites must equal 


N * 1 x 10 18 . 


In the second case, where the only effect of the optical isomers is to 
introduce asymmetry into the environment of the "heavy" atom, the expression 
in equation (6) now equals \ 


AE « 2N6e 


A 

site* 


AE a 2N x 10" 23 au. 


from which we get 
N “ 1 x io 21 . 
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If we assume en area of about 10 A 2 per site, which le a reasonable value for 
a representative SiOi* unit in criatobalite, for example, then the minimum 
area that the single crystal surface must equal fcr a one-percent differential 
adsorption is given by 

A min a * !0 3 c ® 2 (^ rst case) , 

A . “lx 10 6 cm 2 (second case) . 
min 

It is clear from equation (11) that a collection of smaller crystals whose 
total area equals A m<n would not yield the one-percent differentiation since 
the percentage differentiation depends on the number of sites per crystal 
and not on the number of crystals. 

IV. Conclusion 

A cooperative effect is possible on a continuous surface with many 
available adsorption sites. Thus, even though the weak-field Interaction 
terms are almost vanishingly small for one center, a cooperative effect can 
greatly magnify these miniscule energies. Results from a naive statistical 
mechanical treatment show that this amplification feature indeed can lead to 
the resolution of enantiomers by achiral surfaces. For purposes of laboratory 
demonstration, however, the minimum area that an achiral surface must have 
in order to effect a one-percent differentiation is sufficiently large, from 
estimates made in section III, that any experiment would be extremely dif- 
ficult to execute properly. Since, as shown in the text, the quantity of 
discontinuous crystal surfaces does not enhance the percentage of differentia- 
tion, the standard procedure of using large quantities of small particles is 
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not a viable technique for obtaining a detectable (>1%) differentiation of 
optical isomers on an achiral surface. This latter observation may explain, 
to some degree, the ambiguous and often conflicting results obtained to date 
when such experiments have been attempted. 

i 
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